Wave packet dynamics in chains with delayed electronic nonlinear response 
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We study the dynamics of one electron wave packet in a chain with a non-adiabatic electron- 
phonon interaction. The electron-phonon couphng is taken into account in the time-dependent 
Schrodinger equation by a delayed cubic nonlinearity. In the limit of an adiabatic coupling, the 
self-trapping phenomenon occurs when the nonlinearity parameter exceeds a critical value of the 
order of the band width. We show that a weaker nonlinearity is required to produce self-trapping 
in the regime of short delay times. However, this trend is reversed for slow nonlinear responses, 
resulting in a reentrant phase-diagram. In slowly responding media, self-trapping only takes place 
for very strong nonlinearities. 
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The study of the physical mechanisms involved in 
transport phenomena taking place in nonlinear systems 
is a fundamental issue in solid state physics. Concern- 
ing the electronic transport, nonlinearity arises from the 
interaction between electrons and lattice vibrations [l], 
0, S 0, S S 0]- this context, the discrete nonhn- 
ear Schrodinger equation (DNLSE) effectively describes 
the influence of lattice vibrations on the electron dynam- 
ics. The most important property associated with the 
DNLSE is the self-trapping which occurs when the non- 
linearity parameter exceeds a critical value of the order of 
the band width [J, [3, 0|- In this regime, an initially local- 
ized electronic wave packet does not spread continuously 
over the lattice. Therefore, the probability of finding the 
electron at its initial site remains finite in the long-time 
limit. 

In low-dimensional systems, the effect of nonlinear- 
ity seems to be dominant over the role played by dis- 
order [1, 0, [13, III- Recently, the spreading of an ini- 
tially localized wave packet in two nonlinear chains with 
disorder was studied Considering a discrete nonlin- 
ear Schrodinger and quartic Klein-Gordon equations with 
disorder, it was proved that the second moment and the 
participation number of the wave packet do not diverge 
simultaneously ,8;']. The spreading of an initially local- 
ized wave packet in a one-dimensional discrete nonlinear 
Schrodinger lattice with disorder was also recently stud- 
ied [9] . It was observed that the Anderson localization is 
destroyed and a subdiffusive dynamics takes place above 
a certain critical nonlinearity strength. Moreover, an- 
alytical and numerical calculations for a reduced Fermi- 
Pasta-Ulam chain demonstrated that energy localization 
does not require more than one conserved quantity [lo| . 

From the experimental point of view, the interplay 
between disorder and nonlinearity was investigated in 
Ref. [11]. The evolution of linear and nonlinear waves 
in coupled optical waveguides patterned on an AlGaAs 
substrate were directly measured. Nonlinear perturba- 
tions enhance localization of linear waves while induce 



sion to exponential localization was observed. Within a 
more general scope, the study of wave propagation in 
nonlinear systems is a very interesting problem with sev- 
eral applications in many branches of physics. For ex- 
ample, by using a model with diagonal and off-diagonal 
electron-phonon-like terms, the role played by the degree 
of integrability and the resulting extreme event statistic 
were recently studied |l2i] . The main motivation was to 
explain the sudden appearance sea waves of extremely 
large amplitude in relatively calm seas. 

Usually, the nonlinearity produced by the electron- 
phonon interaction is assumed to be instantaneous. How- 
ever, this effective nonlinearity is limited by the finite re- 
sponse time of the medium and corrections are expected 
in the non-adiabatic regime. Actually, the relaxation of 
the nonlinearity is known to have a deep influence on 
the electronic wave packet dynamics [ij, [3, 15, [3]. 
Using non-adiabatic nonlinear models, V.M. Kenkre et 
al 13, 13] considered the problem of a two-level system 
driven by a relaxation process of the nonlinearity. It 
was shown that there is a stationary self-trapping regime 
and a coexistence of static and dynamical transitions for 
certain degrees of nonlinearity and relaxation time [l3| . 
These features contrast with the persistent oscillatory dy- 
namics displayed by the adiabatic nonlinear dimer model 
in both delocalized and self-trapped regimes 15|. More- 
over, it was shown that localized stationary states of 
a dimer model are destroyed above a critical tempera- 
ture 



14| 



delocalization of the nonlinear ones [ll| . In the presence 
of disorder, a transition from ballistic wave packet expan- 



In this letter, we address to the question regarding 
the influence of non-instantaneous nonlinearity on the 
long distance electronic transport in linear chains. This 
physical scenario is particularly suited to the study of 
electronic transport in macromolecules, such as DNA, on 
which the electron-phonon coupling is expected to play 
a relevant role [17|]. Here, we will effectively take into 
account the non-adiabatic nature of the electron-phonon 
interaction by introducing a delayed third-order nonlin- 
earity in the one-electron time-dependent Schrodinger 
equation. Delayed differential equations usually describe 
physical systems presenting a finite response time jl8[. 
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Within a tight-binding approach, the time-evolution of 
the coefficients of the wave vector expanded in the local- 
ized orbitals basis = X]j takes the form 

ic,{t) = Vic,+i{t) + c,_i(t)) - x\cjit - T)pc,(0, (1) 

where we used units of ?i = 1 and considered the on- 
site potentials = 0. V is the hopping integral be- 
tween nearest-neighbor sites, and x is a nonlinear param- 
eter which is proportional to the local electron-phonon 
coupling ^16,] . The non-adiabatic nature of the electron- 
phonon interaction is represented by a delayed contribu- 
tion to the on-site potential proportional to \cj{t — t)^ 
where t is the typical response time. To analyze the wave 
packet propagation, we solve Eq.([T]) using the fourth- 
order Runge-Kutta method to obtain the temporal evo- 
lution of an initially localized wave packet {\'^{t = 0)) — 
SjCj(O)lj) with Cj{0) = (^jjo)- We follow the time evo- 
lution of the amplitude of the wave function at the initial 
site, calculating the so called return probability (19,] . 



Ro{t) = \cM' 



(2) 



Usually, the electron escapes from its initial position 
when the amplitude Cjg(t) vanishes as t evolves. Con- 
versely, the amplitude remains finite for a self-trapped 
wave packet. In order to have a better description of 
the wave packet dynamics, we also compute the time- 
dependent participation function defined as 



Pit) = i/J2\cM 



(3) 



The participation function P{t) gives an estimate of the 
number of sites over which the wave packet is spread at 
time t. In the long-time regime, its scaling behavior can 
also be used to distinguish between localized and delo- 
calized wave packets [l9|. In particular, the asymptotic 
participation number becomes size-independent for local- 
ized wave packets. On the other hand, it scales linearly 
with the chain size in the delocalized regime. 

In our numerical simulations, we set the parameter 
V as the energy unit. Further, we considered the open 
boundary condition with the wave packet initially located 
at the chain center. However, the main following analy- 
sis regarding the self-trapped/delocalized transition and 
phase diagram holds for the periodic boundary condition 
as well. We assume cj {t <0) = for any site j. This con- 
dition corresponds to the situation on which the electron 
is injected in the system at time i = 0, with the elec- 
tronic density being null for negative times. In fig. [IJa), 
we report the time-evolution of the return probability on 
a chain with N = 2 x 10'^ sites for the case of instanta- 
neous nonlinearity (r = 0) for distinct strengths of the 
nonlinear parameter x- The plots show the usual be- 
havior exhibited by the DNLSE with adiabatic electron- 
phonon interaction [3]. For strong nonlinearities x > 3.5, 
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FIG. 1: (a) Return probability Ro{t) for x = 3 and X — 5 
in a chain with N — 2 x 10^ sites with an instantaneous 
nonlinearity. Ro{t) asymptotes a constant value for x > 3.5, 
corresponding to a self-trapped wave packet. It continuously 
decay for x < 3.5 as the wave packet becomes delocalized. 
(b) Ro{t) for X — 5 with different values of t. Even for such 
strong nonlinearity, Ro{t ^ oo) becomes of the order of 
for T = 1. This breakdown of self-trapping is a delay induced 
transition, (c) The participation number P{t) for r = 0.7 
and r = 1 for x = 5. A long delay time allows for the spread 
of an initially localized wave packet. The inset shows the 
linear finite-size scaling of the asymptotic participation num- 
ber for r = 1, typical of delocalized states. (d)(Color online) 
Snapshots of the wave packet in the initial transient regime 
{t = 10), within the plateau (t = 80, 150) on which the width 
of the wave packet around the peaks is roughly constant in 
time, and after the plateau {t — 250) showing the smearing 
out of the wave- front. 



the return probability i?o(0 approaches to a constant 
value in the a long-time limit, thus indicating that the 
wave packet is localized P, 0]- This behavior is the so 
called self-trapping, on which the electron remains local- 
ized around its initial position due to its coupling with 
the lattice vibrations. The critical nonlinearity is of the 
order of the energy bandwidth because x accounts for 
the effective on-site potential felt by the initially local- 
ized electron. If x is much larger than the energy band- 
width, this state can not be efficiently mixed with the 
other energy states. For x < 3.5, the asymptotic return 
probability is vanishingly small (of the order of 1/iV), 
which characterizes the regime of extended states. In 
Fig. [ijb), we focus on the effect of the nonlinearity de- 
lay time T on the self-trapping transition. We consider 
an electron-phonon coupling sufficiently strong to induce 
self-trapping in the adiabatic limit. We observe that the 
gradual increase in the delay time leads to a decrease of 
the return probability. This delocalization process, cor- 
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responding to a breakdown of self-trapping, is a delay- 
induced transition. In Fig. [TJc), we show data for the 
time dependent participation number P{t) in this strong 
nonlinear regime (x = 5). The self-trapping phenomenon 
is indeed destroyed when the delay time is increased, 
with the asymptotic participation number becoming of 
the order of the chain size. It is interesting to notice 
that the participation number exhibits an intermediate 
plateau in the delocalized regime. The time evolution 
of P{t) in this regime reflects the own wave packet dy- 
namics [see Fig. [TJd)]. During the initial stage, the wave 
packet splits from its single peak structure to a structure 
with two peaks that move in opposite direction 20] . The 
participation number grows almost linearly during this 
initial process. After this transient, the wave-fronts dis- 
play a soliton-likc behavior, with the width around the 
peaks being roughly time-independent [ij . In this regime, 
the participation number becomes time-independent be- 
cause the main contribution for it comes from the region 
around the wave packet peaks. Further, there is a dy- 
namical transition to a regime on which the wave-front 
peaks lose their soliton-like behavior. The width of the 
wave-front start to increase in time and, consequently, 
the participation number revivals its growth. It is in- 
teresting to stress that a dynamical transition has also 
been found in nonlinear dimers with non-adiabatic non- 



linearity [13| . Finally, the participation number achieves 
its final saturation when the wave-fronts reach the chain 
boundaries. The linear finite-size scaling, shown in the 
inset, clearly characterizes the delocalized nature of the 
asymptotic wave packet in this regime. 

Data for the long-time return probability Ro{t — > oo) 
versus the delay time t are shown in Fig. [2] for nonlinear- 
ity strengths ranging from x = 2.5 up to 8. The reported 
values are averages in the statistically stationary regime 
of the return probability. These values are finite and 
independent of the chain size when self-trapping takes 
place. On the other hand, the asymptotic return prob- 
ability vanishes as in the delocalized regime. For 
X < 2.5, the asymptotic return probability approaches to 
zero irrespective to the delay time. In contrast with the 
delocalization transition induced by large delay times, a 
short delay can promote self-trapping [Ro{t ~* oo) > 0] 
for nonlinearities slightly below the adiabatic critical one 
2.5 < X < 3.5. The physical reason for the decreas- 
ing of the critical nonlinear strength in the regime of 
fast nonlinear responses is that the nonlinearity comes 
from the electron coupling with the electronic density at 
a previous time. With the electron still relatively local- 
ized around its initial position in the beginning of the 
wave packet time-evolution, the retarded electronic den- 
sity at the initial site is larger than the instantaneous 
one. Therefore, a smaller nonlinear strength is required 
to promote the self-trapping in comparison with the in- 
stantaneous case. Therefore, our results indicate that 
the self-trapping regime has a reentrant character tak- 
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FIG. 2: Ro{t —> oo) versus the delay time r with x = 2.5 
up to 8. A breakdown of the self-trapping phenomenon 
takes place for large delay times. For intermediate electron- 
phonon interactions 2.5 < x < 3.5, we observe that a short 
delay can induce self-trapping [Ro{t (x>) > 0] which is 
absent in the limit of instantaneous nonlinearity. In this 
range of nonlinear couplings, a sequence of delocalized— >self- 
trapped— >delocalized regimes is depicted as the delay time is 
increased. 




FIG. 3: Ro{t ^ oo) versus the strength of the electron-phonon 
coupling X for r = up to 0.8. A delayed nonlinearity can 
induce self-trapping even for nonlinear couplings below the 
critical value in the adiabatic limit x = 3.5. 



ing place in a finite range of short delay times. In 
Fig. [3] we plot some numerical results for Ro{t oo) 
versus the strength of the electron-phonon coupling x- 
Here, we considered delay times ranging from r = up 
to 0.8. These results give further support to the above 
described picture: (i) for large delays the self-trapping 
phenomenon disappears (see results for r = 0.8) in the 
range of nonlinearity strengths shown; (ii) for small but 
finite delays, self-trapping starts below the critical non- 
linearity of the adiabatic limit. It is interesting to no- 
tice that the return probability reaches a maximum at 
a finite nonlinear coupling in the regime of short delay 
times. This non-monotonic behavior also contrasts with 
the continuous increase of the return probability for the 
case of an instantaneous nonlinearity. Finally, we plot 
the phase diagram in the x ^ plane in Fig. [H The 
reentrant character of the transition is reflected by the 
non-monotonic dependence of Xc as the delay time is in- 
creased. The extended states were characterized using 
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FIG. 4: Phase diagram in the X ^ ''' parameter space show- 
ing the delocahzed/self-trapped transition induced by a time 
delayed electron-phonon interaction. The regime of delo- 
cahzed wave packets was characterized using the criterion 
Ra{t oo) ^ l/N and P(t ^ oo) ix N . The critical non- 
linearity has a non-monotonic behavior in the regime of short 
response times. In the inset, we show the critical line in the 
regime of longer response times at which much stronger non- 
linearities are required to promote self-trapping. 

the criterion i?o(i oo) « and P{t — > oo) oc N. 
The slight decrease of the critical nonlinear strength for 
short delay times is replaced by a fast increase at longer 
delays, as shown in the inset. Actually, the wave packet 
spreads substantially before the nonlinearity takes ac- 
tion in slowly responding media. As a consequence, a 
strong nonlinearity is required to stop the subsequent 
wave packet spread. 

In summary, we considered the problem of one electron 
moving on a chain under the influence of a delayed third- 
order nonlinearity which effectively incorporates the ef- 
fect a non-adiabatic electron-phonon interaction. In the 
case of an instantaneous nonlinear response, this system 
presents a delocalized/self-trapped transition at a crit- 
ical nonlinear strength x — 3-5 P, We found that 
for the case of long delay times, self-trapping takes place 
only for much stronger nonlinearities. As a consequence, 
the electronic transport retains its metallic character in 
slowly responding media for a wide range of nonlinear 
strengths. On the other hand, we report a reentrant be- 
havior of the self-trapping transition point in the regime 
of fast nonlinear responses with self-trapping occurring 
even below the critical nonlinear strength of instanta- 
neously responding media. In a finite range of nonlinear- 
ities below the adiabatic critical value, a short delay time 
can promote the stabilization of self-trapping. A second 
transition takes place at a longer delay time above which 
the wave packet becomes again delocalized. 

It is important to stress that the presently reported 
delay-induced self-trapping/delocalization transition is of 
a nature distinct from the one taking place in dimer- 
like structures While in the later the delay pro- 

motes the decay of the oscillations in the population 
levels and a transition from two distinct oscillatory dy- 



namical regimes, the presently reported transition deals 
with the long-distance electron transport. The here re- 
ported phenomenon shall be inherent to dynamical sys- 
tems described by delayed nonlinear differential equa- 
tions. Therefore the influence of the nonlinear response 
time on the electronic wave packet dynamics is expected 
to have a similar counterpart in general (acoustic, electro- 
magnetic) wave transport in nonlinear media. It would 
be interesting to have further studies along this direction. 
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